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Abstract 

This paper reports the initial steps in the devel¬ 
opment of WaterLab, an ambitious experimental 
facility for the testing of new cyber-physical tech¬ 
nologies in drinking water distribution networks 
(DWDN). WaterLab’s initial focus is on wireless 
control networks and on data-based, distributed 
anomaly detection over wireless sensor networks. 
The former can be used to control the hydraulic 
properties of a DWDN, while the latter can be used 
for in-situ detection and isolation of contamination 
and hydraulic faults. 

1 Introduction 

Drinking water production and distribution are 
tightly regulated activities. In most countries, wa¬ 
ter utilities must guarantee minimum levels of wa¬ 
ter quality by keeping the concentration of partic¬ 
ular water contaminants (e.g., bacteria 0 ) below 
levels harmful to consumers [23] . Water quality is 
traditionally monitored by analyzing water samples 
in laboratories every few days or weeks. Although 
there is growing interest in supplementing this rel¬ 
atively slow process with data acquired from new, 
real-time water quality sensor systems [2 EH, the 
testing of this and other new cyber-physical tech¬ 
nologies (e.g., micro-turbines for energy recovery 
m or data-driven anomaly detection 0) in work¬ 
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ing drinking water distribution networks (DWDNs) 
is severely limited due to safety regulations. 

This paper reports the initial steps in the de¬ 
velopment of WaterLab, an experimental facility 
specifically designed to address this limitation. Wa¬ 
terLab will be equipped with a flexible, reconfig- 
urable hydraulic network that, although smaller 
than real DWDNs, will provide the same hydraulic 
conditions (in terms of time constants, flow pat¬ 
terns, etc.) found in practice. This will be 
achieved via an advanced hydraulics control sys¬ 
tem deployed over a state-of-the-art wireless sen¬ 
sor/actuator network composed of hydraulic sen¬ 
sors, motorized valves, and ad-hoc wireless commu¬ 
nication nodes. WaterLab is expected to help test¬ 
ing several new cyber-physical systems technologies 
starting with wireless control networks (WCNs), a 
control paradigm new to the water community, and 
data-based, distributed anomaly detection algo¬ 
rithms executed by wireless communication nodes. 
The facility’s main elements and our initial analysis 
are described in Section 2 next. Section 3 contains 
our conclusions. 


2 WaterLab’s Main Elements 

2.1 Hydraulic Network 

As the simplified diagram in Figure [l] shows, Wa- 
terLab’s hydraulic network is composed of a supply 
loop (green, bottom), a DWDN (top, blue), and 
a set of experimental hydraulic devices (red, left). 
The supply loop is designed to provide an almost 
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Figure 1: WaterLab’s main features. 


constant water pressure differential to the DWDN 
by establishing a high water flow rate across a 
“pressure” valve. The DWDN mimics a residen¬ 
tial water network. The water flow in each of its 
branches will be independently adjusted to match 
the time-varying flow patterns found in typical res¬ 
idential DWDNs. In addition, the DWDN will be 
equipped with manual valves to allow testing un¬ 
der different DWDN topologies. Finally, the ex¬ 
perimental hydraulic devices are add-ons that will 
allow for, among others, the testing of anomaly de¬ 
tection algorithms for pipe-breaks and contamina¬ 
tion faults. The whole setup is expected to be at 
least 20 meters long and to be stacked vertically 
over ground. 

An initial analysis of water consumption data 
from a water utility’s DWDN revealed a daily wa¬ 
ter consumption pattern similar to that sketched 
in Figure [2] As this figure shows, the highest wa¬ 
ter consumption occurs during the morning and 
the evening (when people shower, prepare break¬ 
fast, prepare dinner, etc.), and the lowest con¬ 
sumption occurs at night (when most people are 
asleep). Clearly, there is considerable daily fluctu¬ 
ation. The data also revealed that near the con¬ 
sumption points the most frequent pipe diameters 
are between 40mm and 50mm, and that the flow 
presents slow time constants and is mainly lami¬ 
nar. The setup design and modeling (see below) 
are based on these conclusions. 


2.1.1 Flow Modelling 

Due to the afore-described flow conditions, it was 
concluded that the water inertial effects, the fluid 
compressibility and the pipe’s elastic behaviour 
could be disregarded. As a consequence, the 


flow in the pipes can be accurately modeled by 
one-dimensional rigid water-column equations m , 
which are derived from the Navier-Stokes equa¬ 
tions after certain simplifications (see uni) and are 
extensively used by researchers and practitioners 
EH EDI IB]. According to this approach, the flow in 
every pipe equipped with a valve (see Figure [TJ can 
be modeled by 
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where q (m 3 /s) denotes the pipe flow, g the acceler¬ 
ation of gravity (m/s 2 ), A (m 2 ) is the cross section 
of the pipe, H u (m) is the upstream pressure head, 
H d (m) is the downstream pressure head, L (m) is 
the pipe’s length , / is the friction factor, D (m) 
denotes the pipe’s diameter, and k v is the valve’s 
head loss coefficient. In this equation, the first term 
in the right hand side (RHS) denotes the change in 
flow due to the pressure gradient (assumed to be 
constant as explained before). The second RHS 
term is the change in flow due to friction (which is 
generally a nonliner function of q ). The last RHS 
term denotes the pressure loss due to the valve. In 
this term, k v is a time-variable coefficient that de¬ 
pends on the opening position of the (controllable) 
valve. Thus, k v is the control action variable. In 
the laminar regime / oc 1/g, so the equation above 
can be simplified as follows: 
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where v is the kinematic viscosity (m 2 / s). Finally, 
if the flow is unidirectional (a valid assumption 
when there is only one source of water pressure, as 
is the case in WaterLab), the rigid water-column 
equations further simplify to 


q = JC — Cq — q 2 u , (1) 

where /C = gA Hu ~ Hd > 0. C = > 0 and u = 

2 ^;. This equation is well suited for the derivation 
of hydraulic control laws as explained below (see 
also [6]). 


2.2 Wireless Sensor/Actuator Net¬ 
works 
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WaterLab will be equipped with (among others) 
flow rate, pressure, temperature, conductivity, and 
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Figure 2: Sketch of a typical daily water consump¬ 
tion pattern (red dashed line). The benefits of 
the proposed water flow control law are also shown 
(black dot-dashed line). 

turbidity sensors (“S” in Figure [l]) and motorized 
valves (“V” in Figure [l]). All sensors and valves 
will be interconnected through a wireless commu¬ 
nication network based on an ad-hoc platform (e.g., 
Telosb [26]) that can establish two-way communi¬ 
cation using a variety of communication protocols 
(e.g., ZigBee, gossiping, etc.). The network could 
then be used to transport sensor data to a cen¬ 
tralized location as is done in traditional DWDN 
decision support systems. Alternatively, the wire¬ 
less nodes themselves could be used to process the 
data and implement decentralized algorithms for 
automatic control or anomaly detection purposes 
as explained next. 

2.2.1 Wireless Control Network 

A wireless control network (WCN) is a distributed 
control algorithm executed by the communication 
nodes of a wireless sensor/actuator network IS. 
WaterLab will use a WCN to establish desired wa¬ 
ter flow patterns in its DWDN (see Section 2.1). 
The WCN design requires two steps: First, a con¬ 
trol law is derived based on the rigid water column 
equations using systems and control theory. Sec¬ 
ond, the control law is parallelized. That is, its 
computation is divided into smaller parts that are 
then distributed among the wireless nodes. Our ini¬ 
tial analysis shows that good performance can be 
obtained by using Sontag’s Formula [15]. This con¬ 
trol technique can take into account the valves’ lim¬ 
itations (e.g., saturation), attenuates disturbances, 


and can be easily distributed among the nodes in 
the form of, for instance, a neural network. A sum¬ 
mary of this approach is given next. 

2.2.2 Flow Control via Sontag’s formula 

Consider again the typical daily water consumption 
pattern sketched in Figure [2j Since it is periodic, 
it can be well approximated by a truncated Fourier 
series of the form: 
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R(t) = £ bi cos (iujt) + di sin(icjt), t > 0, (2) 

2=1 

where the constants uj > 0, and G M, 

i = 1, 2, 3, are such that R(t) > 0 for every t > 0. 
The purpose of the control system is to force the 
flow in the pipes to follow the daily flow pattern 
shown in Figure [2] by changing the valve positions 
appropriately. This idea is made more precise next. 

Control Law 

Let e = q — R denote the error between the actual 
and the desired flow in a pipe. It follows from 0 
that the error dynamics are given by 

e = q — R 

= JC — Cq — q 2 u — R 
= JC-£(e + R)-R-(e + R) 2 u 
= /C — £*(c “l - R ) — R — (r 1 R) 2 u — (o • R^ 2 v , 
: = + g(t,e)v (3) 

where v = u — u, for some average valve opening 
position u. The control law v is a real function 
of time, £, and of e (i.e., v : M + xR g R) and 
should be such that e(t) —X 0 as t —X oc for 
all initial error values e(0). Controls laws that 
satisfy the aforementioned conditions are called 
continuous stabilizing control laws and are (in 
general) difficult to implement in practice. A less 
conservative approach, is to design a control law 
such that \e(t)\ —X e as t —X oc for an arbitrarily 
small e > 0, provided that |e(0)| < S, where S can 
be chosen to be arbitrarily large. That is, provided 
that the initial error is bounded, such contin¬ 
uous practical stabilizing control laws can drive 
the error arbitrarily close to zero over time (see 0)- 
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Sontag’s Formula 


Consider the control Lyapunov function m 

V : M + xl4 M + 

(£, e) i —y e 2 /2, 


and let L/V = and = ^g(t,e). It can 
be shown that 


v = 


L g V 

0 


,L g V(t,e) ± 0 
, otherwise. 

(4) 


is a continuous stabilizing feedback control law 
for system ©> D3- Since this control law is 
continuous, the closed-loop system can attenuate 
disturbances such as small measurement noise, 
small actuator errors, small external disturbances, 
and small modeling errors (see [22] Lemma 1.2], 
also 0 Q3] and references therein). In addition, 
note that due to the Weierstrass Approximation 
Theorem, © can be arbitrarily approximated 
using (sufficiently high order) polynomials over e 
and t. Moreover, it can be shown that such an 
approximated controller is a continuous practical 
stabilizing control law. 


Example 

Figure [2] shows the result of applying the flow 
controller to a 100m long pipe under a (constant) 
head differential of 1.85m, assuming a pipe diam¬ 
eter of 46mm and water viscosity of 10 -6 m 2 /s. 
Note that the controller is able to quickly steer the 
almost zero initial flow (8.6 x 10 -6 m 3 /s) to the 
value required by the reference signal R(t). 


2.2.3 Data-Based, Distributed Anomaly 
Detection 

The wireless nodes can also be used to detect un¬ 
expected behaviors or events, such as leaks or in¬ 
creased contamination, commonly referred to as 
anomalies [5]. There is an extensive body of re¬ 
search on this subject, and several anomaly detec¬ 
tion methods are already available to recognize un¬ 
expected events in sensor data. Anomaly detec¬ 
tion is commonly performed at centralized storage 


and processing facilities. It entails creating models 
of the sensor data through either mathematical or 
statistical mean, and then using these models to 
forecast the behavior of future data points. Dis¬ 
crepancies between the new data and their fore¬ 
casts indicate the presence of possible anomalies. 
Examples of these methods can be found in data- 
mining systems mm, data center management 
[25], security card access [1], or network intrusion 
detection m- 

The use of centralized method has several draw¬ 
backs. For instance, in a wireless sensor net¬ 
work, all measured data would need to be trans¬ 
mitted to a selected centralized location. Hence, 
there is growing interest in distributing anomaly 
detection methods in embedded network nodes. 
This approach can provide advantages in terms of 
both data reduction (e.g., only possible anomalous 
data points are transmitted) and resource efficiency 
(e.g., lower bandwidth communication channels). 

WaterLab will provide the opportunity to test 
several new distributed anomaly detection algo¬ 
rithms based on machine learning (ML) techniques. 
Particular attention will be given to online recur¬ 
sive regression techniques that are able to fit both 
linear [4] and non-linear [3] models to sensor data 
using limited-resource embedded systems. These 
techniques do not require supervision (i.e., they run 
autonomously) and are well suited to process data 
from environments or processes for which little a 
priori knowledge is available (e.g., the difference 
between “normal” and “anomalous” behavior). In 
such cases, techniques based on learning or self¬ 
adaptation are easier to apply than their model- 
driven counterparts. 

Since its operating conditions will be well-known 
and modelled, WaterLab will provide the opportu¬ 
nity to compare the relative benefits of both data- 
based and model-based anomaly detection meth¬ 
ods. Moreover, WaterLab will also be used to test 
for the first time in the DWDN context anomaly 
detection methods (currently under development) 
based on locally-combined sensor data from multi¬ 
ple sensors nodes. 


3 Conclusions and Next Steps 

This paper summarizes the main features of Wa¬ 
terLab, an ambitious test facility composed of a 
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scaled-down but dynamically equivalent hydraulic 
network together with a wireless sensor/actuator 
network (WSAN). WaterLab is designed to care¬ 
fully mimic the working conditions of real drinking 
water distribution networks (DWDNs), thus allow¬ 
ing the testing of new cyber-physical technologies 
that could otherwise not be tested due to safety 
regulations in commercial DWDNs. 

WaterLab initial analysis phase is now com¬ 
pleted. The next steps of the project include the 
careful recreation of WaterLab’s hydraulic network 
(together with its control system) in a real-time 
simulation environment (e.g., Simulink) and the 
construction of a prototype hydraulic network (in¬ 
cluding a supply loop, a simple DWDN, and a small 
WSAN). 
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